Dynamics of fibers in a wide microchannel 
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Dynamics of single flexible non-Brownian fibers, tumbling in a Poiseuille flow between two parallel 
solid plane walls, is studied with the use of the hydromultipole numerical code, based on the 
multipole expansion of the Stokes equations, corrected for lubrication. It is shown that for a wide 
range of the system parameters, the migration rate towards the middle plane of the channel increases 
for fibers, which are closer to a wall, or are more fiexible (less stiff), or are longer. The faster motion 
towards the channel center is accompanied by a slower translation along the flow and a larger fiber 
deformation. 



I. INTRODUCTION 

Flexible fibers dynamics has been the subject of in- 
terest in a myriad of problems ranging from industrial 
applications around the paper making process or water 

Sirification down to the level of single DNA molecules 
There are articles, which concentrate on the be- 
havior of suspensions and interaction between fibers lead- 
ing to grouping and entanglement and the flocculation 
phenomena 0, Others look closer at physical mecha- 
nisms of lateral, cross-stream migration, fiber dynamics 
and stretching in the fiow, due to shear or wall effects 

The single fiber problem has attracted a lot of atten- 
tion in recent years due to relevance to many new ex- 
perimental results conducted in microfluidic devices and 
implications on DNA dynamics and its biological role on 
gene regulation [ll|. DNA flexibility is extremely im- 
portant for its functioning. Looping of the fibers allows 
proteins bounded at distant sites to act synergistically 

. An insight into these important functions is now pro- 
vided by experimental techniques which have evolved to 
a point enabling one to directly visualize and manipulate 
the conformational changes irnparted by a fiuid fiow on 
single, flexible DNA polymers [l2l, fl3| . Such experiments 
have shown the conformational dynamics of individual, 
fiexible polymers in steady shear usually exhibit tempo- 
ral fiuctuations described as end-over-end tumbling of the 
molecule [l|. 

Many studies have also concentrated on the phe- 
nomenon of cross-stream migration in confined flows due 
to its relevance in nano- and microfluidic devices and 
fundamental practical significance [l3 - [l9j . Lateral mi- 
gration, which in low-Reynolds-number dynamics is ab- 
sent for rigid particles, is found in the evolution of fiowing 
fiexible fibers due to dynamical changes of shape while 
translating along the channel. It is accepted that confine- 
ment has a crucial impact on the behavior of single poly- 
mers [13, HO] by effects of spatially varying shear rate and 
hydrodynamic interactions with the impenetrable walls. 
In general, the flow properties of fiexible fibers have been 
shown to be highly influenced by hydrodynamic interac- 



tion, both within the fiber and with the system bound- 
aries [2l| . 

Although single-particle experiments have provided a 
leap in the visualization of dynamics of fibers in flows, 
in particular the dynamics of a DNA molecule or a 
nanofiber(2^. theoretical concepts and numerical mod- 
els are often still the key to understanding the under- 
lying physical mechanisms. The so called bead model 
|23| has now been accepted to reproduce the hydrody- 
namic properties of a fiexible fiber. Each pair of bonded 
spheres can stretch and bend, by changing bond distance 
and bond angle, respectively. The strength of bonding, 
or fiexibility of the fiber model, is defined typically by two 
parameters of stretching and bending. Altering these pa- 
rameters, the property of the fiber model can be changed 
from fully rigid to fiexible. It is used widely to model 
fiexible fibers in various simulation techniques ran ging 
from Monte Carlo simulations [IJ] to Brownian 2^ 2^ 
and molecular dynamics (27j . 

In this paper, we describe a single fiber with the use 
of the bead model, and the hydromultipole numeri- 
cal code 1281, w hich implements the theoretical multipole 
method |29l . |30[ of calculating hydrodynamic interactions 
between the particles in Stokes flows. We assume that 
the Brownian motion of the fiber segments is negligible, 
and approximate the channel boundaries by two infinite 
parallel solid walls. We evaluate hydrodynamic interac- 
tions between the segments of a fiexible fiber, and apply 
the single- wall superposition ^31] to take into account the 
influence of the channel boundaries, following the theo- 
retical approach and the numerical procedure outlined in 
Ref. (2^. 

Here we focus on the guiding question how to increase 
the migration rate of the fiber towards the middle plane 
of the channel, by an adequate choice of the system pa- 
rameters. In Sec. ini we describe the system and explain 
our model. In Sec. lIIIl we present the results for the fiber 
motion and its shape evolution. The main conclusions are 
outlined in Sec. IIVI 
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II. MOTION OF FIBERS IN POISEUILLE 
FLOW: THEORETICAL MODEL 

A. Problem 

In the paper, we consider a single mobile and flexible 
microfibcr immersed in a fluid flow. We assume that ve- 
locity of the flow is very small, with the Reynolds number 
Re <C 1, small enough to neglect inertia effects. Here Re 
is the product of a particle velocity and its size, divided 
by the fluid kinematic viscosity. We assume that the 
Peclet number is large, Pe ^ 1, and Brownian motions 
are irrelevant. Here Pe is the product of a fiber veloc- 
ity and its size, divided by the diffusion constant. We 
assume also the absence of unsteady effects. For such 
a system, the hydrodynamic friction forces and torques 
exerted by the moving fibers on the fluid are equal to the 
external forces and torques imposed on them. 

For the system specified above, the fluid velocityv and 
pressure p satisfy the stationary Stokes equations j32l . [33j . 

TyV^v- Vp = 0, (1) 
V-v - 0, (2) 

where rj is the fluid dynamic shear viscosity. 

The fluid is confined inside a microchannel made of 
two parallel infinite solid walls at z = and z = h, as 
shown on Fig. [T] The walls are non-deformable. 




FIG. 1: The system. 

The external fluid flow vq inside the channel is a 
Poiseuille flow, 

vo^4:V„,z{h- z)/h'^x. (3) 

Note that Vm is the the maximum fluid velocity, attained 
at z = h/2. The above set of partial differential equations 
is supplemented by the stick boundary conditions at the 
surface of the flber and at the hard walls, which confine 
the fluid. At inflnity, we assume that v — Vq = 0. 

In this paper, the fiber is initially straight and aligned 
with the fiow, and we detect its motion and shape defor- 
mation. To this goal, we first need to specify the elastic 



and bending forces which keep together all the flber seg- 
ments. These are analyzed in the succeeding subsection. 

B. Model of a fiber: elastic and bending forces 

To determine the fiber motion, the bead model is 
used [1^ . Each fiber strand is constructed out of N solid 
non-deformable spherical particles of diameter d which 
can move with respect to each other. The relative mo- 
tion of the beads is constrained by, and results from elas- 
tic and bending forces, which are discretized as described 
in Ref. 

The elastic (extension) energy between all the neigh- 
boring beads is given by the Hook's law, 

1 ^ 2 

E''^^~kJ2{k-lo) , (4) 

1=2 

where Iq is the equilibrium distance between beads, k is 
the Hooke's constant and 

k = \Ul (5) 

where 

*i = Ti-n_i, (6) 
and Vi is the position of the bead i center (see Fig. [5]). 




FIG. 2: Bead-spring model of a fiber [3^ . Each bead feels 
non-hydrodynamic forces due to the stretching and bending. 

The extensional force exerted on the bead i is then 
equal minus the gradient of the elastic energy, 

F[ = -V,E- = -k{k - + Hh+i - (7) 

where Vi is the derivative with respect to and ti = 
U/k. 

The bending energy in case of a continuous model of a 
fiber is given by the integral, Q 

where A is the bending stiffness. 
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The discretized version of Eq. 
model [33, is 



|), describing the bead 



where 



A ^ 

^ ^ /, (1 - t;+i ■ u) , 

^ 2 — 1 

(9) 



1 for 2<i<N 
for l,iV. 



(10) 



The corresponding bending force = —'ViE^ is then 
equal to 



h-\- 



h 

fi+l . 
U+1 



■t, 



A. 

h+i 



—i- ■ i 



H+1 • ti+2 



fi+l , 
U+1 



H+2 I 



(11) 



Summarizing, the total non-hydrodynamic force ex- 
erted on each bead i — 1, ...,N consists of the extensional 
force and the bending force, 

F = F[ + F,K (12) 

Note that the total non-hydrodynamic force applied to 
all the beads of the fiber is zero, 



JV 



(13) 



From now on we use dimensionless variables. Distances 
are normalized by the bead diameter d, velocities by the 
maximal velocity w,„ of the ambient Poiseuille flow, and 
the forces by 



Time is given in units of 

We introduce dimensionless parameters: 
k = kd/fo, 
lo — h/d, 
A = A/ifod'). 



(14) 

(15) 

(16) 
(17) 
(18) 



The parameter A describes how large is the bending en- 
ergy (or bending force) of the fiber, if the external fluid 
flow ([3)) is fixed. A larger value of A means a stiffer 
(less flexible) fiber, and this is why A is called a bend- 
ing stiffness parameter. At A = the fiber may bend 
freely without the use of energy. In fact, A is the ratio 
of the bending force to a hydrodynamic force related to 
the ambient fiow amplitude v„i ■ In a similar way, k is the 
ratio of the elastic force to a hydrodynamic force related 
to the ambient fiow amplitude Vm- 



Fiber dynamics: theoretical and numerical 
methods 



Assume first that the fiber is fixed at a given configura- 
tion, and that the forces Fq = (Fqi, Fqat) and torques 
To — (Tqi, Tqat) exerted by the fiow on each bead 
are known. For a mobile fiber, Fo and Tq are used to 
determine the translational and rotational velocities of 
all the fiber beads, i — l,...,N, U — (Ui,...,Ujv) and 
fl = (ill, rjAr), respectively. U and depend linearly 
on the elastic and bending forces, F = (Fi, ...,F]y), with 
Fj given by Eq. ([T^ . and on Fq and Tq, 



F + Fo 
To 



(19) 



where /x is the mobility matrix [3(/| . The mobility matrix 
depends on the instantaneous positions r = (ri,...,rjv) 
of all fiber segments. For a given configuration r, values 
of Fq, To and fi are determined by the HYDROMUL- 
TIPOLE numerical algorithm [2^, with the wall effects 
evaluated by the single- wall superposition [2^, [sH . The 
evolution of the fiber is determined by numerically solv- 
ing the set of coupled differential equations. 



dr/dt = U, 



(20) 



where U is given by Eq. ([T9|) . An adaptive Runge-Kutta 
method, implemented in the FORTRAN code is used. 



III. RESULTS 

A. Basic information 

In our study, fibers are initially aligned with the fiow. 
In our frame of reference (see Fig. [1]), the initial con- 
figuration Vi = {xi,yi,Zi) of all fiber segments satisfies 
yi = 0,Zi = zo,i = 1,...,A^. This ensures that at all 
times t, the fiber stays in the plane y = 0. In the per- 
formed numerical simulations we fix the fiowing values: 



k = 10, 
lo = 1.05, 
h = 50. 



(21) 
(22) 
(23) 



A large value of k ensures that the fiber holds an almost 
constant length of 



(TV -1)^0 + 1, 



(24) 



with the maximal relative increase from the equilibrium 
value ((24| of 7% or less, and the minimal relative decrease 
not exceeding 4%. 

In this work, we analyze the dependence of the fiber 
dynamics on its initial position zq across the channel, its 
bending stiffness A and the number of beads N. The goal 
is to determine the range of the parameters which corre- 
sponds to the fastest fiber migration toward the middle 
plane of the channel. 
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B. Typical evolution of a single fiber 

The motion of a single fiber strand in Poiseuille flow 
shows a generic behavior. Initially aligned with the flow, 
it slowly begins a tumbling motion, which is almost pe- 
riodic. Almost, not exactly, because the fiber has a ten- 
dency to migrate across the flow. Fig. [3] shows subse- 
quent snapshots from a typical fiber dynamics. 



To quantify evolution of the fiber position and shape, 
we define the following quantities: 

• r is the tumbling time, defined as the time between 
the two consecutive crossing of the fiber end-to-end 
vector (which links the centers of the first and the 
last beads) with the plane perpendicular to the flow 
direction. The tumbling time depends on A, N and 
the fiber position z across the channel. 

• Xrnid(t) and Zmid{t) dcuotc the X and z position 
components of the middle segment of the fiber, if 
N is odd, and the average position of the beads 
with labels N/2 and N/2 -I- 1, if iV is even. 

• a;Af(^)and ZM{t) denote the x and z components 
of the fiber center-of-mass position. The center is 
calculated as the simple arithmetic mean over po- 
sitions of all fiber segments. 

• 5x{t) and 5z{t) denote the x oi z components of 
the difference between the positions of the last and 
the first segments of the fiber. 



z 

25 — 




Dependence on initial position across the 
channel 



The dependence of fiber dynamics on the initial posi- 
tion zq across the channel has been analyzed for a fiber 
with = 20 segments and the stiffness A = 0.5. 

In Fig. m the position Zmid of the middle point of 
the fiber is shown as a function of time, for different 
initial positions zq. For zq ~ 5,10,15 the evolution of 



N = 20 A = 0.5 
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e 10 
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FIG. 4: Position z^id of the middle section of the fiber as 
function of time, for different initial positions zo of the fiber 
aligned with the flow. 



FIG. 3: Typical evolution of the fiber shape. Number of beads 
= 20, the initial fiber position zo — 15 and the bending 
stiffness A = 0.5. The tumbling time between the first and 
the second flipping moment t = 2585. 



fibers is characterized by tumbling and migration towards 
the center of the channel. Both phenomena depend on 
the instantaneous distance from the wall Zmid(t). Fibers 
placed exactly on the line of symmetry of the channel 
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N = 20 A = 0.5 
lOOOi . . . 




FIG. 5: Difference between the fiber center-of-mass position 
XM{t) along the flow and the reference position of the fluid 
element x ■ vo{zM{t))t, moving with the Poiseuille flow, for 
different initial positions zq of the fiber. 



move along vi^ith the flow without changing their relative 
configuration. 

Fig. [5] shows the difference between the center of mass 
position XMit) of the fiber along the flow, and the posi- 
tion of an element of the fluid moving with the ambient 
Poiseuille flow calculated at ZM{t)- The fiber center-of- 
mass lags behind the ambient flow For the considered 
initial positions, this effect is larger for a smaller value 
of zo, i.e. for a fiber which migrates faster towards the 
middle plane of the channel. 

In Fig. [6l 5z is plotted versus 5x for different initial 
positions zq of the fiber. Each trajectory corresponds 
to two tumbling times: it starts from the 2nd flipping 
moment and ends on the 4th flipping moment of the fiber. 
As expected, 5x\sz=a practically does not depend on zq - 
length of fibers aligned with the flow is insensitive to their 
position across the channel. However, 5z\sx=q is larger 



N = 20 A = 0.5 
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FIG. 6: Difference 5z in position of the last and first section 
of the fiber in the perpendicular to the fiow [z) direction as 
a function of the difference 5x in position of the last and the 
first sections of the fiber in the flow (x) direction, for different 
initial positions zq. 



for larger distances from the wall, where fibers deform 
less. For a smaller distance from the wall, the fiber shape 
is more deformed in comparison to the initially straight 
line. Rigid non-deformable rods would just rotate, and 
their trajectories would be circles. 

Conclusions from this subsection are the following. 

1. For z > 5, fibers perform a tumbling motion. The 
tumbling time increases when the distance from the 
wall for z becomes larger. 

2. Migration rate of fibers towards the center of the 
channel is larger for a smaller distance z from the 
wall (at least for z > 5). 

3. Flipping tends to hinder the motion of the fiber 
along the flow, in comparison to the ambient 
Poiseuille flow. 

4. The faster fiber migration across the channel, the 
larger hindrance of its motion along the channel. 

5. For a larger distance from the wall, the shape of a 
tumbling fiber is less deformed. 

D. Dependence on length (number of segments) 

The equilibrium length of a fiber is given in terms of 
the number of beads N used to model it, and the dis- 
tance Zq between them, see Eq. In this section, we 
analyze how the fiber dynamics depends on the number 
of beads N . We fix the stiffness parameter A = 0.5 and 
the initial position zq = 10- 

Fig. [7] shows the position Zmid as a function of time, 
for fibers made of a different number of segments TV. All 
fibers are initially placed at z = 10 and aligned with 
the flow. For a sufficiently large number of beads TV, 
we observe that longer fibers migrate faster toward the 
middle of the channel, and their tumbling time is slightly 
longer. The shortest fiber, with N — 10, tumbles and 
migrates slowly towards the wall. 

Fig. Oshows the difference between the fiber center-of- 
mass position along the flow, and the reference position 
of the fluid element moving with the ambient Poiseuille 
flow velocity, see Eq. ([3]), calculated at ZM{t)- Longer 
fibers move slower than the ambient flow, and the short- 
est fiber with iV = 10 moves faster. The fibers which 
migrate towards (away from) the wall are faster (slower) 
than the ambient Poiseuille fiow. In agreement with the 
observations from the previous subsection, hindrance of 
the motion along the fiow is correlated with the enhance- 
ment of the migration rate towards the middle plane of 
the channel. 

Fig. [5] shows 5x as a function of time, for fibers made 
of different numbers of beads N . This plot clearly shows 
the lengthening of the tumbling time with the increasing 
length of the fiber. Note that for a shorter fiber, the 
times it takes to fiip is much shorter, compared to the 
time it spends stretched out in the fiow. In particular. 
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A=0.5 Zg=10 



A = 0.5Zjj = 10 




X 10 



FIG. 7: Position Zmid of the middle section of the fiber as a 
function of time for different numbers A*' of segments in the 
fiber aUgned initially with the flow at zq = 10. 




X 10 



FIG. 8: Difference between the fiber center-of-mass position 
XM{t) along the flow, and the reference position x-vo{zM{t))t 
of the fluid, for fibers made of different numbers of beads A'^. 



for N = 10, the plot of 6x{t) consists of almost horizontal 
sections connected by almost vertical segments. 

Fig. [TU] shows 5z as a function of 6x, for different 
numbers of segments in the fiber. The short fiber with 
N = 10 segments does not fold significantly, and the 
plotted evolution curve is almost a circle. The longer the 
fiber is, the more it bends during tumbling. Note that 
Sz\sx=a is smaller for longer fibers - the fiber ends are 
closer to each other at tumbling moments. 

Conclusions from this subsection are the following. 

1. Shorter (longer) fibers move along the flow faster 
(slower) than the ambient flow, and migrate to- 
wards (away from) the wall. 

2. Longer fibers migrate faster towards the middle of 
the channel and move slower along the flow. 




-50, 



X 10 



FIG. 9: Difference Sx in position of the last and the first 
section of the fiber in the flow {x) direction as a function of 
time, for fibers made of different numbers of beads A'^. Initially 
the fiber is aligned with the fiow at zq = 10. 



A = 0.5Zjj = 10 



10 

-10 



N=10 N=15 N= 


20 - 


— N=30 


1 ; 
















f- 1 : 




























.1 ■ 1 ■ I ■ ■ 1 ■ I ■ I 1 


1 I ■ I ■ 1 I ■ I ■ 1 ■ I ■ I 1 ■ I ■ I 1 





-30 



-20 



-10 




5x 



10 



20 



30 



FIG. 10: Difference 5z in position of the last and first section 
of the fiber in the perpendicular to the flow (z) direction as 
a function of the difference Sx in position of the last and 
the first sections of the fiber in the fiow (x) direction, for 
different numbers A*" of particles in the fiber. Initially the 
fiber is aligned with the flow at zq = 10. 



3. Longer fibers take much larger fraction of the tum- 
bling time to flip, and bend over much more during 
that process. 

4. Tumbling time increases with the length of the 
flber. 



E. Dependence on bending stiffness 

The dependence of fiber dynamics on the bending stiff- 
ness A has been analyzed for a fiber with = 20 seg- 
ments and the initial position zq — 5 across the channel. 

Fig. [TT] shows the position Zmid as a function of time 
for different values of the bending stiffness A. The fibers 
move with tumbling towards the central plane of the 
channel. For the bending stiffness A = 0.25,0.5 and 1, 
migration rate of these fibers across the channel increases 
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N = 20 z„ 



N = 20 = 5 




X 10 



FIG. 11: Position of tiie middle section of the fiber as a 
function of time, for different values of the bending stiffness A. 



when the stiffness A of the fiber becomes smaller. How- 
ever, the migration rate ioi A = 2 is larger than for ^=1. 

Concurrently, as shown in Fig 1121 the fibers which mi- 
grate towards the middle plane of the channel, are slower 
than the ambient Poiseuille flow. In agreement with the 
previous observations, more pronounced hindrance of the 
motion along the flow is correlated with the enhancement 
of the migration rate towards the middle plane of the 
channel. 

N = 20 z„ = 5 




FIG. 12: Difference between the center of mass position xai (t) 
of the fiber along the flow and the reference position of the 
fluid element x-V(){zM{t))t, for different values of the bending 
stiffness A. 

Fig. [T51 shows 5z as a function of Sx for different flexi- 
bilities of fibers at = 20 and initial position at zq — 5. 
As expected, Sz\sx=o is larger for more stiff fibers (which 
deform less), while Sx\sz=o does not depend on A (length 
of fibers aligned with the flow is insensitive to A). 

Conclusions from this subsection are the following. 

1. Tumbling fibers migrating towards the middle 
plane of the channel lag behind the ambient 
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FIG. 13: Difference 5z in position of the last and first section 
of the fiber in the perpendicular to the flow (z) direction as 
a function of the difference Sx in position of the last and the 
first sections of the fiber in the flow (x) direction, for different 
values of the bending stiffness A. 



Poiseuille flow. 

2. The faster fiber migration toward the middle plane 
of the channel, the larger hindrance of its motion 
along the channel. 

3. Larger values of A (greater stiffness) correspond to 
slightly larger tumbling times. 

4. The higher the bending stiffness A, the less folded 
is the fiber when tumbling. 



Dependence on the magnitude of the Poiseuille 
flow 



In this section, we analyze the dependence of fiber dy- 
namics on the magnitude of the Poiseuille flow. For the 
comparison, we have chosen N = 20 and zq = 5. If we 
decrease the amplitude Vm of the Poiseuille flow dSj by 
a factor a, then it follows that the relative fiber stiff- 
ness A and the relative elastic parameter k (which are 
inversely proportional to Vm) increase by a factor of a, 
see Eqs. (HH) and 

Here we take as the reference A — 0.25 and k = 5. 
The decrease of the flow magnitude a — 2,4,8 times, 
corresponds to {A,k) = (0.5, 10), (1,20), (2,40), respec- 
tively. In this section, we discuss the resulting change of 
the fiber dynamics. 

Fig. [14] shows Zmid as a function of time, initial po- 
sition of the fiber and the flow intensity. The num- 
ber of segments = 20. Comparing the curves for 
iA,k) = (0.25,5), (0.5,10), (1,20), we see that the 
smaller flow intensity causes slower migration towards 
the middle of the channel. However, the migration rate 
is again enhanced for still smaller flows corresponding to 
{A, k) = (2, 40). As shown in Fig[T5l the fibers which mi- 
grate towards the middle plane of the channel, are slower 
than the ambient Poiseuille flow. In agreement with the 
previous observations, more pronounced hindrance of the 
motion along the flow is correlated with the enhancement 
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N = 20 z„ = 5 




0.5 1 1.5 2 2.5 3 



' x10' 

FIG. 14: Position of the fiber middle-point Zmid, as a function 
of time, for decreasing flow amplitudes, i.e. increasing values 
of {A, k). 

of the migration rate towards the middle plane of the 
chamiel. 

N = 20 z„ = 5 



1000 




FIG. 15: Difference between the center of mass position XM{t) 
of the fiber along the flow and the reference position of the 
fluid element x ■vo{zAf{t))t, as a function of time, for decreas- 
ing flow amplitudes, i.e. increasing values of {A, k). 

An increase of the flow amplitude (i.e. decrease of A 
and k) causes also a decrease of the value of 6z\sx=o, and 
more generally, a larger fiber deformation, as shown in 
Fig. [H 



N = 20 = 5 

— A=0.25 k=5 — A=0.5 k=10 — A=1 k=20 — A=2 k=40 
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FIG. 16: Difference Sz in position of the last and first section 
of the fiber in the perpendicular to the flow (z) direction as 
functions of the difference Sx in position of the last and first 
section of the fiber in the flow (x) direction. 

Conclusions from this subsection are the following. 

1. Increased flow amplitude leads to a larger fiber de- 
formation; fibers tend to keep folded along the flow 
rather than to stretch across the flow. 

2. The faster fiber migration toward the middle plane 
of the channel, the larger hindrance of its motion 
along the channel. 

IV. CONCLUSIONS 

Dynamics of single flexible fibers entrained by the 
Poiseuille fiow in a wide channel made of two parallel 
solid plane walls has been evaluated with the use of ac- 
curate numerical code hydromultipole. The fibers are 
initially aligned with the flow, at a distance zq from the 
wall. For moderate values of zq, and of the fiber bending 
stiffness A and its aspect ratio N, the following generic 
behavior has been observed. 

In our computations, fiber deformation increases when 
they are closer to a wall, or consist of a larger number of 
beads, or are more flexible, or are entrained by a faster 
ambient flow. Moreover, faster migration rate towards 
the middle plane of the channel is accompanied by a 
slower translation along the flow, with velocities smaller 
than the ambient Poiseuille flow. 

For a wide (but not the whole) range of values for 
zq, a and N, migration rate towards the channel middle 
plane increases for fibers, which are closer to a wall, or 
more easy to bend (less stiff), or longer, or entrained by 
a faster ambient flow. A more detailed studies are under 
progress, and will be published elsewhere. 
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